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ABSTRACT 


An arbitrarily orientable Yagi-Uda arrav antenna was 
modeled, and a computer simulation πη to obtain the input 
impedance, gain pattern and front-to-back ratio of various 
arrays. The model made provisions for the antenna to be 
operated over either a lossy ground plane or aboard a ship 
in seas of specified state. Quick solution turn-around, 
with CRT display, enabled relatively rapid optimization of 
numerous arrays. 

Theory, resultant optimal designs and performances, 


photographs, and program listing are included. 
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I. INTRODUCTION 


A. BACKGROUND OF THE STUDY 

The most concise equation that describes the radiated 
electric field of an antenna assumes that the antenna 
Operates either in free space or over a perfectly conducting 
ground plane. Equations in this category are simple ones 
that have been typically used in text books, and voluminous 
experimental data have been obtained through their use. 
Emphasis has been upon simplification of equations so that 
problems could be solved manually with a minimum of rigor. 
Equations such as these are admittedly inaccurate because 
they omit the component of radiation produced by ground 
reflection or if not omitted the ground surface is considered 
to be a perfect reflector.” 

In reality the simplest antenna operates over a lossy 
ground and established at least a direct and a ground 
reflected wave. When simple free space equations are modi- 


fied for two path propagation over a lossy grounc plane and 


lrhe ground wave is composed of a space component and 
a surface component. A particularly important point mav be 
made that for the horizontally polarized wave in the VHF/UHF 
range the strength of the space wave greatly exceeds that of 
the surface wave so as to render the surface wave negligible. 
This then allows an accurate field strength to be calculated 
using equations which assume only a two-path propagation 
model--the direct wave and the lossy ground-reflected wave. 
Henry R. Reed and Carl M. Russell, UHF Propagation, (New 


York: ohn Wiley & Sons, Inc, , 1953), p. 172: 
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solutions are obtained by computer with a Cathode Ray Tube 
(CRT) display of the gain patterns, interesting results 


accrue. 


BS STATEMENT OF THE PROBLEM 

The initial problem is to obtain and verify the equations 
required to describe the electric field of a multi-element 
linear array, particularly of the Yagi-Uda type. A recent 
work which provided graphic computer solutions to single 
element antennas is the basis of this study, and as such 
this study is a follow-on.? Therefore effort centers around 
adapting the existing program to meet the requirements of 


both antenna types. 


en: OBJECTIVES 
The objective of this study is to provide a near real- 

time computer graphic solution of the gain pattern of an 
arbitrarily oriented Yagi-Uda array which is centered above 
two types of planes (with specified e & o): | 

1. the lossy ground plane, and 

2. the lossy ocean plane which rolls and pitches the 

antenna as specified by the sea state. 

Examples of gain patterns and other output parameters are 


illustrated. 


2R. W. Adler and C. B. Robbins, "The Solution and Graphic 
Display of Gain and Patterns for Wire and Linear Antennas in 
the Presence of Lossy Ground", Electrical Engineering Depart- 
ment, Naval Postgraduate School, to be published . 


rn 
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D. SCOPE AND LIMITATIONS 

The following assumptions have been made for the study 
of the Yagi-Uda array: 

l. Propagation is confined to two paths. 

2. The ship-ocean model does not make any provision 
to augment the number of wave paths that turbulent 
seas might produce. 

3. The current on the elements is distributed sinusoi- 

ο νεο ene elements are thin (dr< 1/100). 

4. The elements may be spaced arbitrarily with arbitrary 
lengths. Thickness can be changed. Assignment of 
the element lengths is constrained by the fact that 
if elements are one wavelength then the solutions 
become indeterminant. 

For purposes of testing the resulting equations, various 
array designs were —! The dimensions of the array 
were varied using the method of iterative search to uncover 
optimal horizontal and vertical designs of arrays with and 
without a reflector. 

Optimality is determined by three criteria: 


l. Input Impedance (2. } 


in 
Where it was possible to do so, the array was 
designed to have a reasonably high resistance, e.g. 
> 20% and a reasonably low reactance, e.g. < 10%. 


2. Front to Back Ratio (FBR): 


Within a satisfactory impedance range, the FBR is 
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maximized.? This expression differs in some cases 

from the expression for FBR that is typically used. 
3. Power Gain (G): 

Finally, within the maximum FBR the gain is 

maximized.4 

Photographs are shown of the linear and logarithmic 
results that are obtained from a variety of Yagi-Uda 
designs placed at different heights above ground and ocean 
environments. 

Preliminary tests were made which showed that gain 
varied with change in element thickness, however for the 
results obtained throughout this study the element thick- 
ness to element length ratio remained fixed at 1/200. 


Detailed operating procedures are found in Appendix A. 


JFor purposes of the study the following expression 
for FBR isused 


9 
FBR = G(6,¢) max y! ο $ «180 
G(9,%?)max , 180 <360° 


This differs completely in some cases from the expression 


that is typically used: 
G(9,90) 
O 


, 


FBR = 


“The power gain used is measured with respect to an 
isotropic source. 
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TIE BHBEORY- OF. THESMODEE 


A. ARRAY IMPEDANCE AND ELEMENT CURRENT 

The radiation problem BR res the solution of the 
individual element currents, which will become the basis 
for calculating the gain of the array. The assumptions 
are that the driven element is excited by one volt, and 
that the antenna is thin and therefore the element currents 
are sinusoidal. The equations for the current column 
matrix and for Imax i are 


| EE 
[ty = [2] ἐν] (1) 


li (2) 


max i C D o 
sin B £i 


[2] is the combination of the free space impedance and 


image impedance which is found by 


^ 


[2] - ΟΤΙ (3) 


[2] and [2'] are NxN dimensional--[Z'] being the mutual 
impedance matrix relating the actual and image elements as 
indicated in figure l. It should be noted that the matrices 
are complex. Complex inversion presents no particular 
problem except that the [2] matrix must be arrayed according 


to the basic arithmetic operations that follow: 


14 
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Figure l. Side view of a 3 element 
Yagi-Uda array tilted & 
degrees above the ground. 
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(RJ + 30xD) + (itp) + 0p) = WI WM 


which can be written as two equations 


[R][1g] - [xlL1í] 


[VI re 


(5) 


[x][r.] + [R][I;] imm. 
which can be solved according to equation (1), where 


equation (3) now appears as: 


Ξ [κ] [-x] 
[z] = (6) 
[x] [R] 
[21 is a 2Nx2N matrix. An array of maximum size (5 
elements) requires inversion time on the SDS-9300 computer 
of 5 seconds--roughly one-eighth the total solution time 
needed for the complete aximuth and elevation patterns. 
The reflection factor C in equation (3) is a function 


of the angle which the antenna element makes with the 


ground. C is written as” 


C = e73^' (R, 'cos A' * jR;'sin A'). (7) 


The reflection coefficients Rn ' and RS are the values 


obtained when 90-0 is substituted into equation (8). Theta 


μ. T. Ma and L. C. Walters, Power Gains for Antennas 
over Lossy Plane Ground, ESSA Technical Report. ERL 104-ITS 
74 (U. S. Government Printing Office, Washington, D.C., 
1969). 
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is equal to zero because the coupling which takes place 
between one free-space and the image element (thereby 
producing a mutual impedance Z') occurs with the image 
directly beneath the actual element. Equations (3) and (7) 
show that the horizontally polarized field (^'20^) and the 
vertically polarized field (Δ'-907) give weight to the 
value of [2'] by the values C=R,' and Sab respectively. 
Since the solution to equation (3) is actually independent 
of 0 and $, one solution for this equation satisfies the 
gain expression at any position of observation. 

The values of Rn and Ryr which are also used in the 


gain equation, are a function of the observation angle 9 





as ο. 
β2 
cos 6 - a A 
Rh — le (8a) 
cos 6 + B2 A 
B 
cos 9 - = A 
R; = (8b) 
B 
cos 9 + A 
Ba 
where 
2 3 
ο ο ο (98) 
β2 
ὃν 
=: _— + € 9b 
82 Bl r J acc (2b) 
i πο š 
82 = Ble, - j) ===] | (9) 
FMH 
6rbid. 
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The matrix [Z'] is identical to equation (6) except 
that it is written in terms of [R'] and [x']. The geometric 
orientation of the antenna elements are central to the 
solution for the mutual impedance matrix [Z']. Figure 2 
shows an arbitrarily oriented single element and its 
corresponding image in the ground plane. Although the 
true coordinates of the antenna are in the xyz coordinate 
system, the equation for solving for the impedance fixes 
the primary element along a vertical axis--here it is shown 
to be z''. The secondary element becomes the image element 
and is separated from the primary element by the distances 
Zo and Yo The linkage between the two coordinate systems 


t t ' ` ' ` 
xyz and z y'z' is 0'. The equations for 2g and Yo are 


observed to be 


COS O ZA SN A (10a) 
o aa 2h sin 0 = 2h cos A' (100) 
S,t = S cos 20! = S cos 24' (11a) 
Sy = -S§ sin 28' = -S sin 2A' (llb) 


The solution of self impedance is performed by assuming 
that Z =0 and Yo= 2 r, where the radius r is assumed equal 
ἵο £17200. (This makes 1,/d = 100, which for 2; = 1/4 
satisfies the thin antenna specifications of d < 1/100). 
The mutual impedances of the free-space array, where the 
elements are all parallel, are solved using 2 =0 and Yoryır 
the separation between the two elements in question. Under 
dynamic ship motion, which affects 8' and n, the antenna 


height above the earth plane, the mutual impedance 
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image 
(secondary) 
antenna 


Figure 2. The geometry of an arbitrarily oriented element 
in free space and its image (used for solution 
of mutual impedances Z') 
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matrix [Z'] will be affected according to equations (10) and 
(11), but the self and mutual impedances of the [Z] matrix 


will not be affected since it represents the free-space 


array impedance. 


Impedances are solved according to the following: ! 


Do 


: 1Ὠ(2π in(2 
R = -30 f [ee (Zo+S; ' +Ëp) + =- a (Zo+S, ' P) 
-L3/2 1 2X 2 2) 


-2 sin(27r) l 


' ας δω 
2λ 


22 sin(2rr) SD) . Sin(27rj) . sin(27rj) 


e o^ 
i L | 
sin (2 - Is) 
EE cp) BE EM ds (12a) 
S 
L 
se eos r) : Lp, 
X - -30 —————— (2.152 Se 
τ 
Sum a 
ASA 2 cos(2rr TL 
+ EN (Z +87 '-_B) - S cos ( —P) 
ro 2λ r À 





(25*84') = = x. + paon οο5 (56) 


(Yo tS, ) ia λ 


2 ) | s ns - |s|) 
“οοβ(2πτι) | cos(2mr2) 8. zx Be 


(12b) 





7n. C. Baker, and A. H. Lagrone; "Digital Computation 
of the Mutual Impedance Between Thin Dipoles;" IRE 


Transactions on Antennas and Propagation; March, 1962; 
AP-10, No. 2; pps 172-178. 
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EU qs n* (13a) 
2 . 
r [o° + (2 +S;") 1 (13b) 
αρ 
ri = ie + (2 FS, 455) ] (MSc) 
2 E 
ro = Lo” + (2,+5,'-B1"] (13d) 


Where the distances r, ri, and τ», are the respective 
distances in figure 2 from the bottom, center and top of 
the primary antenna to the differential, ds, on the 


secondary antenna. 


B. THE FAR FIELD 

To solve for the electric field, a point of observation 
is chosen, and each element in the array is viewed. Each 
element contributes to the total field according to its 
vector phase and E ics separation from some 
reference element, the tilt angle of the array, and the 
bearing to the point of observation. Figure 3 shows the 
geometry of a 2 element horizontal Yagi-Uda tilted adegrees. 

The equation of the field produced by one element located 
at the origin is first developed. Then, the presence of 
other elements in the array is taken into account. 

The equation at the point P for an element at the origin 
p jo (t-g) 

I = EM sin[8(2,+x)] e (14) 

Current flowing in the incremental length dx produces a 


field, and there will be a phase associated with it 


¿Í 





Array axis (ZY Plane) 


Y 
ith 
element IN, 
NR 


ο Pi 
y 
FR Z Ce E i 
reflector Jo E 
element 7 ᾷ τα S AM TR 
EA T n, l 
` Si 
ar St Plato 


Figure 3. Yagi-Uda array with horizontal elements and 
a tilted axis, as viewed from a distant point ο 


($':0,0 -90) 
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with respect to point P as the differential length moves 
L | ; : 
from >= to + LS SJ Une devay 1s yx cos v, which is found 


ΕΠΕ dot product of ÓX and OR, as follows: 


OA = ÓX + OR 


where. 

^^ ^ ^ ^ 

OX = yfi sin 6' cos 6' + j sin 6' sin 6' + k cos e'] (15a) 

OR = fisine css +4 sine sing +kcose] ( 15b) 
so that 

OA = x[sine sin 6' cos (¢-9') + cos6 cos Ze (16) 
or 

OA = xcos y (175 
where 

COS jg = sin 0 sin 0' cos (ġ-ọ') + cos8® cos 0' + (18) 


τι is now clear that s = Yr =- x cos V and 


: r-x cos 
io Ss: τ sin[8(2.+x)] Jule (EX SOS S| (19) 
i O; i 
where 
b Sul, X ° O 
Ly i. = (20) 
fee τ.» 


Although there is a phase difference between r and s, at 

points taken far from the origin their magnitudes are 

approximately equal. | 
Observe that by translating the xyz origin as well as 


the off-axis element to the point B located along the 


2:3 





^ 


observation vector r, the magnitude of the new observation 
vector relative to the new origin is BR, which differs in 
length from the old observation vector OR. The amount of 
change equals that which affects the observation vector S. 
Furthermore the angle by definition remains unchanged. 
Eberetore, the expression r —- s = xX cos ¥ holds for every 
radiating element in the array, not just for the element at 


the origin. From this it can be concluded that the integral 


equation for the field is identical for each element in the 


array. 
Following the method used by Kraus:® 
hs | 
ME ο ο ο επι J tax, (21) 
i i 
SA 
Substituting for I and letting B - = r “we see that: 
C 
Es 
j 30 jw (t=) j8xcosV 
EQ; a sin lI, e c sin (82,+8x) e? dx , (22) 
us 1 E 


When this expression is integrated the result is: 


ds 
_ j 60 sin 6 I. ¿Ju (ta 
Es; O A A T [cos (82; cos y) -cos (821) ] 28025) 


E Επι Y 


37 p. Kraus; Antennas; (New York: McGraw-Hill Co. 1950), 
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For an array of elements the field is the superposition of 
the fields of each element taken separately. E, is then 


the sum of Eq; given by 


E = ————————————X 


e EE 60 sin ‘So, [cos (82:cosV) -cos (821) ]. (24) 
$ r sin? ψ 

1=1 
It is recognized that the time term, e"jot has been dropped 
and that f has been substituted for /c. 

The array factor, or the phase delay ÓB introduced by 
the separation between elements, is next taken into account. 
From figure 3, using the reflector as reference, the 
separation between elements as viewed from point P is 


e^ 


OB = OY°OR , where 


OY = jy; cos a + ky; sin α | (25) 
and OR = i sin? COS + j sin 0sin$ + k cose (26) 
gives OB = y; (cosa sin 6 sin ¢ + sinacosé) . (27) 
If cos Q = cos a sin 8 sin ὁ + sina cos 8 (28) 
then OB = y; cos Q (29) 


represents the phase lead of the ith element. 


Having found the slant separation between elements to 
be expressed by equation (29), the current for any element 


in the array may be written as 


I - QJBYicos q _ (30) 


= I : 
O; Max 1 


The phasor sum of currents in the array becomes 


25 


ο. ο. 


m 


E 2 


Se ws NIE 


5 





3 Ἴβγασος 2 
Io = 4 Imax i © . (31) 


ı=1 


Upon substitution of equation (31) into equation (24) the 


result is 





: N 
|j 60e JPrsin $ I ,J8Yicos Q 
WE. - > max i 
- r sin y i=1 
[cos(8?:-cos Y) - cos(B2,)](32) 
where T i is found from equation (2). 


Equation (32) represents the free space radiation of 
the horizontally polarized Yagi-Uda. The equation does not 
consider the ground reflections. Figure 4 shows the direct 
and ground reflected waves, and the geometry shows the dif- 
ference between their path lengths. The element separation, 
written as: 


2h; = 2(h + y, sin a) | (33) 


^ 


is a function of the path difference vector AB, since 
AB = 2h;cos 6. 
Thus 


- = Zeh; cos 
E m Ep + Rh E, e 1 , (34) 


And finally, equations (32) and (34) are combined to produce 


the final space wave in the phi direction: 


3 60 es sin ὁ 


Eum | We 
r sın wv 
-j28h;cos 8 
Tax; I Re 1 ]: 
= . . 
e P Yi€95 fT eos (82 cosy) - cos (82,1) ]. (35) 
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Figure 4. The tilted Yagi-Uda receiving horizontall y 
polarized waves via direct and reflected paths. 


(The array is tilted to present less confusion 
in distinguishing angles) 
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Observe that when cos 2 = sin © sin ç and h;=h equation 
(35) agrees with Ma and Walters equation for the horizontal 
qa ` The equation for the 9 component of the horizon- 
F uy g U ES 15 NOt derived™with the Same rigor since 


it represents the cross polarized radiation, which exists 


only off-axis and is of minor interest. The equation is: 


j 60 e 2*®"cos 8 cos ὁ 


Eg 


E r sin V 
-328h;cos ϐ 
I ΙΕ a ] 
max: V 
ms i 
ο” 24998 “Ecos (82; cos V) - cos(B2¿)]. (36) 


The equation for the vertical dipole is well known. 
If p. is used insteady of R, then equations (31) and (34) 
may be applied to the standard free space equation to arrive 


at the general solution for the vertical array: 


-j8r N 
AS -428h. θ 
h;- dE ΠΝ 
r sin ϐ 1 = T 
j8y;cos Q 
e [cos(82;cos ¥) - cos(82.)] (37) 
Specifically, cos Y = sin f sin 9 and cos € z cos 8 
since % =® = 0 for the vertical Yagi-Uda. 


M.T. πα αμα i r. Wabters, opem hp... 


lOrraus, Joezert, 
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C. THE SHIP-OCEAN MODEL 

Figure 5 shows the effects of the dynamic ship model 
upon the Yagi-Uda. Changes in the sinusoidal waves, 
considering that a roll of eight degrees per sea-state is 
produced and a bow pitch of three degrees per sea-state is 
produced (representative of a light cruiser), will affect 


the angles 9' and 9 as follows: 


e' = 6¢' - (A,cos ¢' + Agsin φ') (38) 
a = ας - AjSin $' + A7COS 4' (39) 
where 
T (wave) sin (course) (40a) 
Ay =  (wave).3 cos (course) (40b) 
and where 
wave = 8(sea state) sin (y) f (40c) 


The effective height of the array will change according 
to 


h = h. cos Aj cos A, (41) 


It should be noticed that the parameter ®' until now 
has been somewhat extraneous for the solution to antennas 
over land, since variation of the observation angle $' and 
the angle $ accomplish the same thing. However the angle 
$' now is significant in specifying the antenna orientation 
with respect to the ship heading, which then determines the 


aspect that the wave presents to the antenna. 





Ships 
y Heading 


> 


\ 





Figure 5. Orientation of the Yagi-Uda 
aboard a ship. Angles A, and Az 
represent roll and bow pitch 
respectively. 
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D. GAIN 


The gain of an antenna is equal tomthe ratio of the 


power intensity to the power density, and is expressed as 


follows: 

G = 
where W. = 
or Wà = 
and W. = 

RT 


The final expression for 


The equation for gain becomes independent of r when the 


4 Wr 


Win 








[1 | R. 
b in 
gain appears as 


r2|Ee,1? + |Ee41* 





(42) 


watts/solid angle 


(43a) 


(43b) 


(44) 


equations (35) and (36) are substituted into the equation 


above. 
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III. RESULTS OF ARRAY PERFORMANCE 


A. RESULTS OVERVIEW 

Basically, four types of arrays were examined for 
optimization: horizontal and vertical with and without 
reflectors. Experiments were conducted to determine the 
design of arrays that afford optimal performance as defined 
in the Scope and Limitations section of the study. Starting 
with the two element array, the parasitic element function- 
ing both as a director and as a reflector, the element 
spacings and lengths were varied to obtain optimization. 
The arrays were lengthened by adding one element at a time, 
up to a total of five elements, each time solving for an 
optimal design. Numerous iterations were required, and the 
time required to optimize an array varied from ten minutes 
for a two element array to slightly over three hours for a 
five element array. (The vertical array was more difficult 
than the horizontal, and the array without reflector was 
more difficult than that with reflector.) Not all of the 
multi-element arrays were optimized according to the rigid 
criteria chosen for this study. Specifically, the vertical 
designs having four and five elements were arrived at by 
relaxation of the impedance threshold criterion thereby 


producing a suboptimal design. 


B. THE TWO ELEMENT ARRAY 
Referring to Table I for the horizontal two-element 


array with a director, maximum FBR is maintainec when, for 
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TABLE I 


OPTIMAL 2-ELEMENT HORIZONTAL ARRAY 
AT VARIOUS ELEMENT SPACINGS USING 
THE PARASITE AS A REFLECTOR AND AS 
A DIRECTOR 


(Gain and Front-to-Back ratio in db) 
(Spacing and length in À) 
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f =146. h=lm. In both cases the length of the driven 
element is 1/2. (D απο ο are director απο σερ ες 


lengths respectively.) (The subscript p denotes peak, 
e.g. the length of the parasitic element was varied until 
the peak FBR was first obtained, then within that value, 
the gain peak was found by making finer adjustments to 
length.) 
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increasing element spacing, there is a corresponding decrease 
in director length. Apparently, best performance occurs when 
spacing is .01\ (.021lm) because this is where the maximum 

FBR appears. However, no mention has been made about the 
affect upon impedance of varying the element spacing. 

In general it was found that the resistance and inductive 
reactance decreased together (from 70 to 40 ohms and from 
40 to 10 ohms respectively) as elements were placed closer 
together. When the spacing closed to .051, the currents in 
the elements evidently became very large, and the value of 
resistance became less than one ohm (too low for convenient 
matching). Thus the FBR and G are insufficient criteria of 
optimality since impedance can fall (in some cases may rise) 
to unreasonable values. Optimality in the case of the di- 

Fe CCOr array in Table I actually occurs at „1X (.205m) 
spacing. 

When the parasitic element is used as a reflector it is 
not necessary to change the reflector element length as the 
spacing between elements is changed in search for the 
maximum FBR. The maximum FBR occurs at .5X (1.03m) for all 
spacings except at .254 (.515m); that is ‚the maximum FBR 
occurred when both element lengths were the same. 

The gain curves shown in figure 6 differ from other 
published experimental and theoretical data (and there are 
fairly wide variations in these results too), inasmuch as 
the value of the peak gain is somewhat higher than that 


typically recorded. This is because published data typically 
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Power Gain over Isotropic db 





X NER 
e~. X. 
Ne FBRN 
Ὃς Χ 
12 ` ^x 
en reflector 
game Tr / 
10 e i `e. . M 
1 "e x 
l Sa | 
| `x director 
| Sof 
8 ,FBR `. 
`e 
\ 
e 
1 
\ 
6 \ 
\ 
\ 
\ 
E 
4 N 
@ 
\ 
\ 
\ 
\ 
\ 
2 M 
` 
\ 
N 
`e 
Element Spacing 
| .Φ .ὦ 


Figure 6. Comparison of gain vs 
element spacing for optimal 2— 
element horizontal arrays. These 

' curvesrepresent datain table l, 
where the array isa reflector — 
driven element combination in one 
case and adirector-drivenelement 
combination inthe other case. 
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does not reflect a lossy ground plane (two-path propagation 
is ignored). Also, the gain of the reflector array exceeds 
that of the director array, which generally disagrees with 
other findings. Again, the explanation for these differences 
is tracable to the assumption that a lossy ground wave is 


being propagated along with the direct wave. 


C. THE MULTI-ELEMENT VERTICAL AND HORIZONTAL ARRAYS 

Tables II and III list the optimal and suboptimal design 
parameters ae the fas °° arravs. Notice that the four and 
five element vertical arrays have relatively wide reflector 
spacings--these are suboptimal designs. By comparison with 
the optimal design that was found for the four element 
Mestalcal abray (where R;, > 209) the FBR is 16.4 db whereas 
the suboptimal design (Rin < 200) produced a FBR of 22.4 db 
-- quite a drop. The two designs differed completely.” * 
Furthermore they responded to slight parameter variations 
differently. 

The four element vertical array was fine-tuned for 
maximum FBR only to discover that the impedance was 
unacceptable; all subsequent fine adjustments, introduced 


for purposes of raising the resistance to a value above 202, 


llc, R. Fry, The Yagi-Uda Aerial--A Short Design Review 
and Bibliography, (Valcartier, Quebec: Canadian Armament 
Research and Development Establishment, May 1966), p. 23. 


l?The optimal design, as compared with the suboptimal 
design shown in Table II, was L;=.49,.47,.45,.43; A222, 
15 απο ΕΠ λ}. 
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caused a drastic reduction in FBR--below 16 db. Thus, this 
particular design was abandoned; another completely different 
design provided an acceptable impedance with a peak fine- 
tuned FBR value of 16.4 db (compare parameters in footnote 

/ 
#12 with the comparable parameters in Table II). 

Except for the two larger vertical arrays (four and five 
elements), all changes that were made for the purpose of 
me casing Ri, or reducingeX,, inevitably improved the FBR 
figure and in none of these cases did the manipulation of 
parameters for purposes of adhering to the impedance 
criterion result in obtaining a lower FBR or G value. 
Because such manipulation did, however, cause a drop in the 
FBR for the two vertical arrays they are considered 
anomalous. 

The optimal horizontal and vertical s have 
tapered elements in the sea e next to the driven element, 
but the last few end-directors of the four and five element 
horizontal array show a slight inverse taper. Writers, 
however, SETA agree that elements and spacings that are 
gradually tapered towards the end of an array will usually 
give best results. 

It is obvious that the optimum parameter measurements 
obtained for the reflector-driver combination (two-element 
array) do not necessarily ensure optimization of successive 


experiments involving the addition of another director. A 


1 δημο optimal arrays include all horizontal plus the two 
and three element vertical arrays. 
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design readjustment of all parameters is required each time 
the array length is changed (e.g. each time an element is 
added or removed). A corollary to this for arrays with more 
than two elements is that the FBR will decrease with the 
addition of a director whose length and spacing are identical 
to those of the preceding director. This was observed in 


testing the uniform array!“ 


at fuyg=10 where the FBR dropped 
from 14.6 db as a four element array to 10 db as a five 
element array. 

With regard to director spacings, perhaps it is more 
than just coincidental that the designs which produced 
optimization of the longer horizontal arrays resulted in 
equispaced director elements. Nothing could be obtained 
to substantiate whether other analysts would agree. 

From Table III, comparison between the horizontal three 
and four element arrays as well as comparison between the 
four and five element arrays shows that the FBR improvement 
is greater than the G improvement as array lengths increase. 
An added advantage, besides increasing FBR and G, is that 
the greater the number of elements used in the optimal array 
the greater the tendency for automatic reactance cancella- 
tion. While the reactive components of the two and three 
element arrays are small they did not become so as a 
coincidence of optimization; manipulation of parameters was 


necessary to obtain reactance decreases in conjunction with 


1Mnesign: L;=.533,.500,.434,.434,.434 and d;=.233,.100, 
.100, @ h=3 (L, d and h are in à). 
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the need to increase resistance (and, as inferred earlier, 
these manipulations in the case of the two and three 


ze 


element arrays produced slight decreases in FBR and G 
Comparison between the optimal: horizontal array and its 
vertical counterpart indicates that the horizontal array 
consistently out-performs the vertical array in FBR, but 
matches it in G. Furthermore, when the optimally designed 
horizontal array is operated in the vertical position and 
vice versa a degradation in performance results. An improve- 
ment of almost 16 db in FBR and 1.5 db in G can be obtained 
when the array is optimized in the position for which it is 
intended to be operated. 
Table IV gives the results of operating the optimal 
array without its reflector; it is interesting to observe 
the importance that the reflector has in determining the 
overall array performance.  Practically no directivity was 
obtained without a reflector even with a four element 


(three director) array. The unusually large G that results 


for the three and four element vertical arrays are anomalies 


6 


l5Fishenden reiterates the popular belief that a 
convenient method of altering the input impedance without 
affecting FBR or G is by varying the spacing between the 
reflector and the driven element from 1/8 to 1/4. Such 
attempts in this study did however noticably affect the 
FBR and G (particularly in the case of the vertical four 
and five element arrays, which were particularly trouble- 
some). R. M. Fishenden and E. R. Wiblin, "Design of Yagi 
Aerials", Proceedings IEE, 96, Pt. 3, (Jan 1949), p. 6. 
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TABLE IV 


THIRTY ARRAYS DESIGNED AS SHOWN IN TABLE II 
BUTSWITUTREFCECTPORS REMOVED 


Horiz o Horrz Vert 8 Vert 
Orientation Orientation 
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and cannot be explained (undoubtedly it is connected with 
the fact that these two arrays are actually suboptimally 
designed to begin with.) | 
Such poor performance is not without remedy, however, 
as the results show in Tables V and VI. Arrays without 
reflectors may also be optimized with considerable improve- 
ment in performance. Unlike the arrays with reflectors, 
the optimal design of arrays without reflectors gives 
identical optimal results regardless of whether the array 
is operated in the horizontal or vertical position.  Inas- 
much as the reactance is much greater than 102, these 
arrays do not represent optimal designs according to the 


purposes of this study, yet none better could be found. 


DIS SELECTED INTERESTING PATTERNS 
Numerous photographs were taken of the following antenna 
patterns (in a few cases comparisons between the linear and 
the log gain plots are shown): 
1. a three element array operating over dry land, fresh 
water and sea water, 
2. athree element array at different heights over land, 
3. several three element arrays operating on a ship 
underway in heavy seas, 
4. ἃ three element array with its axis tilted at various 
angles from the horizon and 
5. some resulting patterns when an array is operated 


off-frequency. 
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TABLE V 
PARAMETERS OF OPTIMAL 30 MHz ARRAYS WITHOUT REFLECTOR 


(Horizontal & Vertical Optimal Arrays Identical) 
(L and d in A) 


2-Element 3-Element 4-Element 


E 
= 
T O 
° 

o 
00 


$ 
O 
Ul 





TABLE VI 
RESULTS OF OPTIMAL ARRAYS WITHOUT REFLECTOR 


(Specifications and parameters same as in Table V) 
(G and FBR in db, R and X in ohms) 
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Figure 7 Optimal 30 MHz 
3-element horizontal array: 
linear plot (a) and log 
plot (b).16 


S. τος. 
FBR, = 22.0 db 
BU - 31.6020 MS In 


Figure 8 Optimal 30 MHz 
4-element horizontal array: 
linear plot (a) and log 
plot (b). 


EL = 1353 ab 
FBR, = EMO ab 
Rin = 36.18 X;, = -2.32 


Figure 9 Optimal 30 MHz 
5-element horizontal array: 
linear plots 12) & (b), and 
109 pilot. (ec): 


= 13.5 db 
27.3 db 


G 
m. 








l6parameters that are not shown in Table II but which _3 
are common to all the 30 MHz arrays are: h+=3, e=5. o=10 


$4,290. 0,285, however figure 9 (c) differs in observation, 
with $-20P and 6-60. 


45 





os 


MN MN 








te) 


Fresh Water Sea Water 


G - 12.7 db 
9.5 db FBRp 
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G 
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Figure 10. Ten MH three element vertical 


array radiating over three types of terrain.!/ 


Figure ll. 150 MH Optimal 
three element horizontal 
array--linear prot (a), and 
log plot (b).1 


G = 12,9 αρ 
FBR, A RUE 
Xin =-13.5 an 





8) : (b) 
Linear Log 





171,¡=.506,.500,.450; d,=.233,.133; he=2.; fmy=10.;e=10., 
o=10"3, dn=90., 0p=84. (land); e=80., 0=2., $5790., θρΞ76. (Η20} ; 
e=80., o= 0-2, Φρ-20., 0p=75. (sea). (L, d and ht are in 2) 


185,2.510,.490,.460; d;-.250,.08; ht-15; e-5.; o-1077; 


$p790.; 05789. (L, d and h£ are in à) 
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a 
10 
G = 8.6 @ =80.0 G 0.7 ο ο ποσα Ὁ 
9.6 P | 
(C) (d) 
haere 140 hr = 2.0 
G = 12.9 @ p= /6.0 G 2813.3 E p^ 93. 
FBR, = 27.0 FBR, = 27.3 


Figure 12. Ten MHz 3-element vertical arrays (a) & (b) 


and horizontal arrays (c) & (d) at different height 
over land.l? (G and FBR in db) 





l?see figure 16 for parameters. (h, is in 4.) 
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ω = 0. ω = 30. w = 90. 


stable pt.roll = 15.4  st.roll - 31.0" 
up pitch = 5.5 dn.pitch = 11.0° 

0' - 71.5 6' = 120.0 

o = 64.0 a = -11.0 

h = 1.82 h = 1.68 


Figure 13. Ten MH three element horizontal array 
aboard ship in rough seas. Sea state = 7 and relative 
Airection of waves = 40. Gain varies. 29 





w = 0. w = 960. 
stable pe Poll = 41.45" 
up.pitch = 7.2" 
0' = 48.5 
g cy 
h - 2.80 h - 2.20 


Figure 14. Thirty MH three element horizontal array 


aboard ship in rough seas. Sea state = 6 and relative 
direction of waves = 60. Gain varies. 





2011-.506,.500,.450; dj=.233,.133; e=80.; o=5.; $=90.; 
6 280.; h.-2. (L, d and ht are in 2) 

21,;,7.51,.50,.46; di-.25,.08; e=80.; o0=5.; $=90.; 
6=85.; Ah¿=3. (L, d and ht are in À) 
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Figure 15. 
ship in rough seas. 


of waves 


40. 


SS 
πι ο ρα ση 


240. 


298 ico dl 
9.6°d.pitch 
ot = 10. 

a = -9.6 
h = 13.2 


150 MHz 3-element horizontal array aboard 





22 
~e=80., 


Sea sta 
Gain varies. 


55 = 6, and relative direction 


Other parameters include L;=.51,.50,.46; d;=.25,.08; 


o=5., 


$=90., 


6=89., 


h=15. 
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96, 0 


12.269 =78.0 G 
= 23.4 FBR 
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(d) 
a = 90.0 a = 90.0 α Ξ 90.0 
G = 3.5@8=76.0 G = 13.3 G = 13.3 
Gp - 7.788 =1.0 FBR = 0.0 FBR = 0.0 
$ = 180.0 $ = 0.0 
9 = 76.0 8 = 50, 


Figure 16. Ten MHz 3-element horizontal array at 
various tilt angles of array axis (bore sight) .23 





2Jvarious antenna parameters are: L;=.506,.500,.450; 
.233,.133; e=10., o=1073, φ-90., 0-276. for (a)-(d), 
=1. 


d; 
h+ 


50 


12.1@9=76.0 G ξ 9.1645-76.0 
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Figure 17. Ten MHz 3-element 
horizontal array operating over 
land at fyy=30. (For antenna 


specifications see the ship-ocean 
model figure 14.) 





fMH=25. Euy=30- Mg ??* 
Figure 18. Thirty MHz 3-element 
horizontal array operating over lan 
above and below frequency for which 


array was designed. (For specifications, 
see the ship-ocean model figure 14.) 





Figure 19. Ten MHz 5-element 
horizontal array over land. 
Comparison of linear(inter) 
plot with log plot. 24 


Gp=14.6 db, FBR5-6.5 db. 


Figure 20. Ten MHz 3-element 
horizontal array over land. 
Comparison of linear (inter) 
plot with log plot. 25 
Gp=12.2 db, ο ο db. 





24various parameters include: L 
2131665; di1—-.25,.167,.167,.167; o -10., 
h-1l. 


34,.500,.466,.466, 
0“; ¢ =90.; 6 =77.; 


¿S5Parameters for this array are: L;=.550,.500,.433; 
d;2.233,.133; all other parameters are the same as in 
footnote 24. 


02 








= 
La κα - ου 


E. SUMMARY 

There is a somewhat higher gain and an increasingly 
higher front-to-back ratio as the element spacing widens 
in the case of the optimal two element horizontal array 
possessing a reflecting parasite as compared to the optimal 
array with a director parasite. 

The vertical multi-element array will produce a high 
front-to-back ratio--with values in the neighborhood of 
those of the horizontal array; however, if resistance is 
to be maintained above 202 it appears that the resulting 
optimal vertical array will not provide as high a front-to- 
back ratio as the optimal horizontal array. The optimal 
design of a multi-element horizontal array provides both 
a maximum front-to-back ratio and a near resonant input 
impedance, whereas the best design of a multi-element 
vertical array rss os sas a compromise between these factors, 
e.g. a trade-off sacrifice exists. 

The optimal design arrived at for a horizontal array 
becomes a suboptimal design when operated vertically, and 
vlce versa. 

Upon adding or eliminating directors from an array which 
has been designed for optimal performance--with the stipula- 
tion that parameters of all directors be kept uniform, the 
performance of the array will no longer be optimal. 

Operation of the horizontal array ever sea water environ- 
ment provides a higher gain with respect to operation over 
land, but over fresh water provides a lower gain than over 


land. 
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Gain is considerably higher in the case of the vertical 
array and slightly higher with the horizontal array when 
they are radiating from a height of 2) as compared with 1). 

Gain varies considerably when the platform upon which 
an antenna operates is in heavy seas. 

The gain of a horizontal array aimed vertically into 
the sky is lower than when it operates with no tiit angle 


whatsoever. 
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IV. RECOMMENDATIONS 


At the beginning of this study the existing program 
which had been written for eight different single element 
antennas permitted the compilation of, at maximum, only 
three different antennas--the computer memory was exceeded 
when more options were added. Now the revised program 
allows all nine antennas to be compiled, with the eight 
unused antenna subroutines residing externally on drum. 
This permits consideration of further expanding the main 
program capabilities inasmuch as the effective computer 
memory has more capacity. 

In regard to expanding the program's capability for 
solving array problems, it is recommended that the graphic 
input portion of the program be extended to allow for more 
element length and space parameters. This would permit 
study of long arrays (arrays having a cumulative Spacing 
that total more than one wavelength). 

It is recommended that an optimization routine be added 
to the main program which would encompass one or more of 
the following criteria: gain, front-to-back ratio, 
resonance of the driven element (255) , thickness of the 
element, and impedance (or standing wave ratio). For this 
purpose, a more powerful computer--which gives faster 


solution time--might be preferable. 
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APPENDIX A 


DESCRIPTION OP PROGRAM OPERATING PROCEDURES 


The program for the system of antennas as originally 
written has been changed to accommodate the Yagi-Uda. 
Accordingly, the program was rewritten in part to obtain 
more effective use of limited information storage in the 
computer.  Redundancies between antenna programs have been 
reduced, the shipboard simulator and the gain pattern 
routine was rewritten to accommodate the Yagi-Uda and to 
save storage, and the graphic text segment was changed to 
satisfy the need for additional parameter inputs and outputs. 
The program remains functional for all antennas listed below. 
Figure 21 shows how the parameters and the output appear to 
the operator at the CRT graphics display console. 

TO operate the program, first ready the Graphics device 
to be used by following simple lab instructions for setting 
up "Gated." The program may be quickly loaded into the 
computer by mounting the program "Dump" tape on tape unit #1 
and loading a few cards--BOOT,AJOB, AAGT, AASSIGN Xl1-MT1A, 
and ARERUN. Next, press the usual operating buttons on the 

9300 console--IDLE, RESET, CLEAR FLAGS, RUN, and CARDS. The 
cards and tape are then read by the computer and within five 
seconds the computer console E print a request for the 
graphics device number that the user desires to interface 
with the computer. The user then responds by typing at the 


IDEV=1* and a carriage return. If a mistake is made in typing 





start over again. If the message was successfully typed but 
for some reason the text did not appear on the graphics device, 
the quick load procedure may be used by pressing the following 
computer buttons: IDLE, RESET, STEP, RUN; then retype the 
message IDEV=1*. The graphics device should now contain the 
program text, providing all equipment is functional. 

Each non-zero parameter entry must consist of a field- 
length of fewer than four numbers, plus a decimal, e.g. 
parameters should not exceed four characters of information. 
Upon completing a parameter entry, press the graphics TTY 
carriage-return and the blinking-light pointer will sequence 
to the location of the next parameter input. If a parameter 
is Zero it need not be typed; press the carriage-return and 
go to the next entry location.” Correction to a typing 
mistake may be made only on the four-character field string 
of the parameter where the light-pointer is located. If the 
error is identified after sequencing to the next location the 
program may either be executed as is if no serious program 
error will result, of if the error can cause abort within the 


SDS-9300, either use the fast-restart?9 procedure or reload 


25In addition to the need for a decimal to be entered 
with all data, a blanked-out space (one that has been skipped) 
is the same as a zero. 


26This procedure is similar to the initial start-up: 
type the "Reset Gated” messages at the graphics teletype, and 
at the SDS-9300 press IDLE, RESET, STEP RUN end in response 
to the console message listed by the SDS-9300 type IDEV=1* 
(device #1). 
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the core "dump" tape again. Elther procedure takes less 


than 30 seconds whereas problem solution-time can take longer. 





Figure 21. Photograph of graphic 
display showing parameter input/output. 





27 this procedure is identical to the initial set-up, in 
that the program dump tape is read into the computer. Press 
the button number "32" at the computer console, and type 
ARERUN at the console. This causes the dump tape to rewind 


and be read in again. Once again the computer should reques 
Idev, and the user should respond by typing IDEV=1*. 


Observe the parameters shown in the photograph. The 
geilowingeis a description of these parameters to assist 


the user in specifying his antenna. 


ANTN: Antenna 
l. Tilted Dipole 
2. Vertical Whip 
3. Vertical Whip with Ground Screen 
4. Inverted L 
5. Sloping Long-wire 
6. Sloping Vee 
7. Horizontal Rhombic 
8. Vertical Half-rhombic 


9. Yagi-Uda 


ELE LENG: Element length (meters) 

If a single element antenna (ANTN 1 through ANTN 8) 
has been specified, the first available parameter space 
should be used for the element length and the other four 
spaces to the right left blank. If a Yagi-Uda (ANTN 9) has 
been chosen, the following order must be observed: 


ELE LENG 


νιν 


Driv 
Directors 


The user may choose to design the array without a reflector, 
in which case the first location at the left must be left 


blank. Any number of directors up to three may be selected, 
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however any blanked directors must be contiguous from the 
right. That is, skipping a director space then entering 

data into the next space, will result in the problem being 
misinterpreted: the first zero director marks the end of 


the array and determines the number of elements present. 


HGT/WDTH: Height and width (meters) 

A single element antenna requires only that the height 
be specified; the last four spaces may be left blank, since 
they will be ignored. If the Yagi-Uda is chosen, select the 
parameters according to the following format: 


HGT/WDTH 


En 2 ae 


Spacings 
The first spacing is for the distance between the reflector 
and driven elements. Subsequent spacings are for the distance 


between respective directors. 


9 
PHIP: φ' Antenna parameter. Normally .¿'=0 will be used 





except under dynamic ship simulation, where the angle will 
specify the array axis relative to the ship's head. See 


figure 5. 





THEP: θ' Antenna parameter. 6 '=0° for a vertical YAGI-UDA 


560 for horizontal. 


FREQ: Frequency (MHz) 
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EPSL: Dielectric constant of earth 





SGMA: Conductivity of earth (mho/m) 


Typical values are: 





i 
t 
1 
I 
! 
t 
4 
: 


UL: Multiplier 

jJ EMubtHbUXres o0 by 1 

2  Multiplies m by .001 
For numbers “< 3 this parameter may be used as a multiplier 
of the diameter of the array elements. The element diameter 
is used in determining the element spacing for computation 
of self impedance. When SMUL » 3 the,expression used for 
Spacing is 

£i 


d | 2 X; : 
—— or 1, since r= i. 
y2 200 200. 


(see the computer program instruction that precedes #1925) 
PHI: 6 Observation parameter. 


MET: 9 Observation parameter. 

Generally the cbserver wishes to position himself at 
the point of maximum gain, and take the azimuth and elevation 
sections there. The first solution of the problem will show 
where this point is located by listing the respective angles 
Under PK) "Peaks. The user may then entergthese values for 


PHI and THET, and solve the problem again. 
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LOG: Log pattern 

T£ left blank linear patterns will result. If "1." 
is entered, a previously entered pattern via light pen 
will be erased. If "2." is entered, a log pattern will 


result. 
ISTH: Store the horizontal pattern 


mov; Store the Vertical pattern 

Leave blank except to store the pattern about to be 
calculated for purposes of future comparison with another 
pattern by means of "double-exposure". Enter "1." to 


store. 


IRCL: Recall a stored pattern 

Left blank until patterns have been stored by a 
previous problem solution using ISTH/ISTV. When comparison 
1s desired with the solution of the problem being entered, 
enter "1". Ensure that if recall is desired ISTH and ISTV 
are left blank, otherwise comparison will be made of the 


Solution with itself. 


DEEN: Array-axis angle 

The tilt angle (or the angle of elevation with the 
ground) of the Yagi-Uda array axis. This is the slope 
angle with the ground for the two sloping antennas (ANTN 5 
and 6) as well as for the vertical half-rhombic (ANTN 8), 
and the minor ET between elements of the horizontal 


phombic (ANTN 7). 
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ISEA: Sea state 

For the tilted dipole, vertical whip, sloping long- 
wire, and Yagi-Uda. Left blank unless iterative solutions 
under dynamic sea conditions are desired. The program is 
written for a sinusoidally varying sea with solutions 
produced at thirty degree intervals. The peak amplitude 
of roll is 8 degrees per sea state, and the peak pitch is 


3 degrees per sea state. 


ICRS: Course of the seas 
Enter the direction of the seas relative to the ship's 


heading. 


Bu. Gain at the observation points specified by ¢, o. 
The peak gain, identified under "PK" is also given. 

The ç and 0 which locate the peak gain are listed adjacent 

to the values unse θη For the problem (discussed above 


under PHI and THET). 


FTBR: Front-to-back ratio 
The ratio of the peak forward-azimuth to peak reverse- 
azimuth gains. (This is not a comparison between the 90 


and 270 degree gain values). See footnote 3. 
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